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Abstract 

Within the generalized high-twist factorization formalism, we express the contribution from multiple parton scat- 
tering and induced gluon radiation to the DY dilepton spectra in terms of nuclear modified effective beam quark 
distribution function. We show that beam quark energy loss is characterized by jet transport parameter q, which is 
related to the local gluon density of the medium. Using the value of q determined from the deeply inelastic scattering 
(DIS) data, we evaluate the nuclear modification factor in the Drell-Yan process in p+A collisions. Effects of parton 
energy loss in the DY spectra are found negligible in the Fermilab experimental data at = 800 GeV relative to 
parton shadowing while the predicted suppression of the DY spectra are significant at E^b = 120 GeV. 
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1. Introduction 

Jet quenching is one of the most important discoveries in heavy-ion collisions at RHIC. It can provide crucial 
information of the properties of the hot and dense matter [1 1, quark gluon plasma (QGP), which is expected to be 
formed under extreme conditions. However, in A+A collisions, both final-state QGP effects, such as jet quenching, 
and initial-state effects in cold nuclear matter can modify the experimental observables. In pursuit of the real signals 
from QGP, one has to study the baseline initial conditions in A+A collisions, which include transverse momentum 
broadening, parton shadowing and parton energy loss in cold nuclei. 

In this study, based on generalized high-twist factorization formalism |2|, we discuss the initial-state effects of par- 
ton energy loss in the Drell-Yan process in p+A collisions due to multiple scattering inside the nuclear medium. The 
calculation is similar to final-state multiple scattering in deeply inelastic scattering (DIS) off large nuclei, where the 
medium modification to the quark fragmentation function caused by multiple scattering and induced gluon bremsstrahlung 
is effectively written as Q, 

b h q {z h ,ri 2 ) = D h ( Zh ,n 2 ) + -^ \ -f -\Ay q ^ g (zJl)D h q (z h /z) + Ay^ gq (zJl)D h g (z h /z)] , (1) 

Zn Jo lj Jzh Z 

where the medium modified splitting functions Ay a ^bc(z, l\) depend on the properties of the medium via the jet 
transport parameter which is related to the local gluon density, 

4n 2 a s C F 

^f(^n) = N2 _ 1 Pa(^n)[xG(x)] xx0 , (2) 
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where is the nuclear density distribution and G(x) is the gluon distribution function in a nucleon. Based on the 
modified fragmentation function from single induced gluon emission, one can derive a medium modified DGLAP 
evolution equation, and determine the value of q from comparisons to the HERMES experimental data [4]. A recent 
study finds a range of the jet transport parameter qo = 0.016 - 0.032GeV 2 /fm in cold nuclei by fitting the HERMES 
data on and K* spectra in SIDIS |5 1 and this value is consistent with the transverse momentum broadening of the 
Drell-Yan dilepton production in p+A collisions |6). In addition, a recent study shows that one can also obtain the 
value of q by studying the nuclear enhancement of transverse momentum imbalance for two back-to-back particle 
production in e+A and p+A collisions J7J- 

2. Parton energy loss in cold nuclei 

In the Drell-Yan process in p+A collisions, the energetic quark coming from the beam hadron will undergo mul- 
tiple scattering with the target remnant before it annihilates with an antiquark from the target nucleus to produce 
dileptons. The contribution from multiple scatterings in nuclear medium can be calculated by the generalized high- 
twist factorization formalism in pQCD [2 1. In such a factorization formalism, one has to expand the hard partonic part 
of the hard interaction in the intrinsic parton transverse momentum. In the expansion, the first term gives rise to the 
normal collinear factorized pQCD results, known as leading twist contributions. The high-order terms in the Taylor 
expansion are known as higher-twist contributions which can be expressed as the convolution of hard partonic parts 
and high-twist matrix elements inside the target nuclei. In general, the twist-4 contribution can be expressed as 



where £ is the momentum fraction carried by the initial parton in the beam hadron and kj the intrinsic transverse 
momentum of the initial gluon in the nuclear target. In our study, we consider two distinctive double scattering 
processes as shown in Figs. 4 and 5 in HI . One is the annihilation-like process with the additional scattering between 
the beam quark and a gluon from the nucleus; the other is the Compton-like process, where the additional scattering 
is between two gluons coming from the beam hadron and target nucleus, respectively. Besides these two double- 
scattering processes, one also needs to consider the interferences between single and triple scatterings. Contributions 
from the first term in the collinear expansion from all the processes mentioned above give the eikonal contributions 
to the single scattering to ensure the gauge invariance of the leading-twist results. The second derivative of the 
hard partonic part with respect to kj gives rise to the twist-4 contributions. However, the contributions from the 
interferences between single and triple scatterings are power suppressed {qj/Q 2 ) as compared to the leading-twist- 
4 contributions from double scattering. Within the generalized high-twist factorization formalism as shown in Eq. 
[3] one can calculate the final leading logarithmic contributions coming from the annihilation-like and Compton-like 
processes. Summing up all the leading logarithmic contributions from single and double scattering, we can express 
the contributions from multiple scattering and induced gluon radiation to the DY dilepton spectra in terms of nuclear 
modified effective beam quark distribution function ||8] , 



with x' is the momentum fraction carried by the beam quark in the projectile hadron. The modified beam quark 
distribution function takes a form very similar to the vacuum bremsstrahlung corrections that lead to the evolution 
equations in pQCD for parton distribution functions. In Eq. (j4j), fq/h(x' ,/j 2 ) is the renormalized twist-two beam quark 
distribution. The nuclear-dependence of the medium-modified beam quark distribution function is implicit through 
the medium-modified splitting functions Ay 9 _, 9g and Ay g _»^ which are defined in Eq. (32) and (37) in [8]. 




(3) 




(4) 
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The medium-modified splitting functions are proportional to twist-4 matrix element T g£j which can be factorized in 
terms of quark and gluon density distribution inside the target nucleus under the assumption that nucleon correlation 
inside nuclei is negligible. The gluon density in turn can be related to quark transport parameter q in the nuclear 
medium, therefore, 



T„Jx, x t ) 



6§q/§/jv(jc) 



na s p A (Q. 



— I d 2 b I df I dy\p A (y\,b)p A (y~ ,b)sm 2 (x,p 

0j_) «J DO xj — DO 



(5) 



where we have used the assumption that the jet transport parameter along the quark jet trajectory is proportional 
to the nuclear density q(y,b) = qop A (y,b)/p A (Q,Q), and we will employ the Woods-Saxon nuclear geometry in our 
numerical calculation. In the definition of T g *, the trigonometric function comes from the LPM interference which 
plays an important role in the small q 2 T region and leads to the quadratic nuclear size dependence of the parton energy 
loss. 

Shown in Fig. [T]are the calculated nuclear modification factors at leading order, 



Ba^ _ Bdo- pA ^i-/dQ 2 dx' _ £ ? J dxf ij / A (x,p 2 )f q / p (x' ,p 2 ,A)H (x,p,q) 
Ac^ ~ AdcT pB ^,-/dQ 2 dx' ~ AZj dxf m (x,p 2 )f q/p (x\p 2 )H (x, p, q) ' 

as a function of x' for the Drell-Yan dilepton production as compared to the Fermilab E866 experimental data J9). 



(6) 
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Figure 1: Ratios of the DY cross section per nucleon in p + A collisions versus x'at E lab = 800 GeV. The shaded bands 
correspond to qa = 0.024 ± 0.008 GeV 2 /fm |5 1. The experimental data are from the Fermilab experiment E866 l9l . 



For fixed invariant mass M, the fractional momentum x = M 2 /x's carried by the target partons becomes smaller 
for large beam quark fractional momentum x' in the kinematic region of the E866 experiment, therefore strong nuclear 
shadowing of the quark distribution inside the target nucleus. This behavior is clearly demonstrated in the comparison 
between our calculation and the experimental data in the kinematic region of the E866 experiment. The dominant 
nuclear modification of the DY cross section is from nuclear shadowing of parton distribution functions inside large 
nuclei as given by the EPS08 [ 10] parameterization. The effect of medium-modified beam quark distribution caused 
by beam quark energy loss leads to further suppression of the DY cross section for large nuclei. However, with the 
quark transport parameter predetermined from the nuclear DIS experiment 0, the suppression due to initial beam 
quark energy loss is quantitatively small. The additional suppression only becomes considerable in large x' region in 
a large nucleus. Since the parameterization of nPDF [ 10 1 from global fitting included DY data, one should include the 
effect of beam parton energy loss in large x' or small x region. 

On the other hand, in p+A collisions at lower beam energy Ei a t, in p + A collisions, the target parton momentum 
fraction x becomes large for moderately large beam parton momentum fraction x', where the effect of shadowing 
should be small. In this kinematic region the fractional parton energy loss will become larger for smaller beam parton 
energy x'Em, 0. Therefore, at fixed DY dilepton mass M and lower beam energy E/ a b, one can disentangle the effect 
of initial-state parton energy loss from the nuclear shadowing. Shown in Fig. [2] are the predictions for the DY cross 
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Figure 2: Predictions for the DY cross section ratios per nucleon in p + A collisions versus x' at E k , b =120 GeV. The 
shaded bands correspond to q = 0.024 ± 0.008 GeV 2 /fm 0. 



section ratios at E/ a i, = 120 GeV in the Fermilab's E906 experiment [11] with invariant dilepton mass M = 4.5 GeV. 
At this lower beam proton energy, the effect of parton shadowing is indeed small as shown by the dashed line in Fig. [2] 
On the other hand, the energy loss effect induced by multiple scattering is significant and the dominant cause for the 
DY suppression shown by the shaded bands in Fig. [2] Therefore, the E906 experiment can provide an unambiguous 
measurement of the effect of initial-state parton energy loss in DY cross section. 

3. Conclusion 

We have discussed the medium modified beam quark distribution function in the DY process in p+A collisions 
due to multiple scattering within the framework of generalized high twist factorization formalism. The medium 
modification to the beam quark distribution depends on the twist-4 matrix element which can be related to jet transport 
parameter q. Using the value of q determined from nuclear DIS data, we evaluate the nuclear modification factor in 
the Drell-Yan process in p+A collisions for different nuclear targets in two kinematic regions corresponding to the 
Fermilab E866 and E906 experiments, respectively. We found that in E866 experiment, the effect of beam parton 
energy loss is only considerable in the large x' region with heavy nuclear targets. However, this effect becomes 
significant for lower beam proton energy as in E906 experiment. 
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